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ABSTRACT 

Middle East respiratory syndrome coronavirus (MERS-CoV) utilizes dipeptidyl peptidase 4 
(DPP4) as an entry receptor. While bat, camel, and human DPP4 support MERS-CoV infection, 
several DPP4 orthologs, including mouse, ferret, hamster, and guinea pig, do not. Previous work 
revealed that glycosylation of mouse DPP4 plays a role in blocking MERS-CoV infection. Here, 
we test whether glycosylation also acts as a determinant of permissivity for ferret, hamster, and 
guinea pig DPP4. We find that while glycosylation plays an important role in these orthologs, 
additional sequence and structural determinants impact their ability to act as functional receptors 
for MERS-CoV. These results provide insight into DPP4 species-specific differences impacting 
MERS-CoV host range and better inform our understanding of virus-receptor interactions 


associated with disease emergence and host susceptibility. 


IMPORTANCE 

MERS-CoV is a recently emerged zoonotic virus that is still circulating in the human population 
with a ~35% mortality rate. With no available vaccines or therapeutics, the study of MERS-CoV 
pathogenesis is crucial its control and prevention. However, in vivo studies are limited because 
MERS-CoV cannot infect wildtype mice due to incompatibilities between the virus spike and the 
mouse host cell receptor, mouse DPP4 (mDPP4). Specifically, mDPP4 has a nonconserved 
glycosylation site that acts as a barrier to MERS-CoV infection. Thus, one mouse model strategy 
has been to modify the mouse genome to remove this glycosylation site. Here, we investigate 
whether glycosylation acts as a barrier to infection for other nonpermissive small animal species, 


namely ferret, guinea pig, and hamster. Understanding the virus-receptor interactions for these 
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DPP4 orthologs will help in the development of additional animal models while also revealing 


species-specific differences impacting MERS-CoV host range. 


INTRODUCTION 

Coronaviruses are a diverse family of viruses that infect a wide range of hosts, including 
both mammalian and avian species. Phylogenetic studies suggest that over the last 800 years, 
several zoonotic coronaviruses have expanded their host range into humans, resulting in four 
antigenically distinct strains that are still circulating in the human population (1-2). Whereas 
most human coronaviruses cause only mild symptoms in healthy adults, two recent emergence 
events have resulted in severe disease in humans. Severe acute respiratory syndrome coronavirus 
(SARS-CoV) emerged from its zoonotic reservoir in 2003 and infected over 8,000 people with a 
9% mortality rate before being controlled by public health measures (3). Although SARS-CoV 
is no longer circulating in the human population, several SARS-like coronaviruses persist in bats, 
can use the human receptor for entry, and are poised for emergence (4-5). In 2012, Middle East 
respiratory syndrome coronavirus (MERS-CoV) emerged from its zoonotic reservoir and 
continues to cause human infection. As of July 2017, MERS-CoV has infected 2,037 people with 
a ~35% mortality rate (6). The exact origins of MERS-CoV remain unclear. However, most data 
suggest that the virus originated in bats and spread into human populations using camels as an 
intermediate host species (7-8). Understanding how MERS-CoV and other coronaviruses evolve 
and spread will allow us to combat MERS-CoV infection while also developing approaches for 
dealing with future pandemic coronavirus strains. 

One of the key determinants of viral host range is the interaction between the virus spike 


protein and the host cell receptor. For MERS-CoV, the host cell receptor is dipeptidyl peptidase 


FJILSVOMGN JO AINN AQ ZL02 ‘82 Ainr uo /Bio"wseiAl//:dyjy wo. pepeojumoq 


(0) 
= 
6 
oO 
(0) 
— 
(72) 
O 
jae 
— 
(OL 
= 
O 
72) 
=) 
(= 
= 
TO 
ve 
[ok 
0) 
O 
O 
~s 


Journal of Virology 


Journal of Virology 


76 


77 


78 


79 


80 


81 


82 


83 


84 


85 


86 


87 


88 


89 


90 


91 


92 


93 


94 


95 


96 


97 


98 


IV (DPP4) (9), a ubiquitously expressed cell surface protein that functions in immune 
homeostasis. Interestingly, whereas MERS-CoV can utilize the bat, camel, and human DPP4 
molecules for entry (10-12), it is unable to infect cells using the DPP4 molecules from traditional 
small animal models, including mice, ferrets, guinea pigs, and hamsters (13-17). The inability of 
MERS-CoV to infect these species in vivo and in vitro is primarily due to spike-receptor 
incompatibilities and not due to other species-specific host cell factors (13, 18-19). Because 
MERS-CoV is unable to infect traditional small animal models, our ability to study MERS-CoV 
pathogenesis is restricted, as is the ability to evaluate drugs and vaccine efficacy. To overcome 
this obstacle, several mouse models have recently been developed in order to study MERS-CoV 
pathogenesis in vivo. These include an adeno-associated virus delivery of human DPP4 (hDPP4) 
(18), overexpression or knock-in of hDPP4 (20-23), and the use of gene editing to generate mice 
carrying a chimeric mouse DPP4 (mDPP4) gene (24). These mouse models greatly improve our 
ability to study MERS-CoV infection, motivating us to investigate the potential for alternate 
species to act as models for MERS-CoV infection. Importantly, the development of vaccines and 
therapeutics relies on positive results from more than one animal species. For MERS-CoV, the 
disparate findings reported for nonhuman primate systems, including the rhesus macaque and 
common marmoset (25-28), emphasize the need for additional animal models. 

In our previous studies, we found that receptor incompatibilities between MERS-CoV 
and mDPP4 can be alleviated with just two amino acid substitutions (13). Introducing the 
mutations A288L and T330R (residues numbered relative to mDPP4) allows mDPP4 to act as a 
functional receptor for MERS-CoV infection. The A288L mutation helps stabilize a hydrophobic 
core in the MERS-CoV receptor binding domain (RBD) shown to be important for proper 


binding to DPP4 (29). Conversely, the T330R mutation knocks out a nonconserved glycosylation 
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site present in mDPP4. This glycosylation site is a substantial barrier to MERS-CoV infection, 
with mDPP4 only acting as a functional receptor when the glycan is absent (19). Taken together, 
these two mutations informed the development of the aforementioned chimeric mDPP4 mouse 
model (24) and improve our understanding of the biochemical mechanisms that govern the 
interactions between DPP4 and the MERS-CoV RBD. 

Here we investigate whether glycosylation acts as a broader determinant of permissivity 
using alternate small animal models including the ferret, hamster, and guinea pig. Overall, we 
find that while glycosylation of these DPP4 orthologs plays an important role in blocking 
MERS-CoV infection, additional species-specific determinants contribute to the inability of each 
of these receptors to support MERS-CoV infection. Further investigation is required to identify 
these determinants and establish whether they lie at the level of the sequence or structure. 
Improving our understanding of the interactions between DPP4 and the MERS-CoV RBD can 
help inform the development of alternate small animal models as well as help reveal how 


coronaviruses are able to emerge into novel species. 


RESULTS 

Susceptibility to MERS-CoV varies among a number of host species, making detailed 
analyses of orthologous DPP4 receptors essential for elucidating fundamental mechanisms that 
govern coronavirus species specificity. DPP4 orthologs are classified as permissive (human, bat, 
and camel) or non-permissive (mouse, ferret, hamster, and guinea pig) (Figure 1A). Based on 
mean fluorescent cell count, hDPP4, bDPP4, and cDPP4 all support higher levels of infection 
compared to nonpermissive orthologs (p < 0.05, Student’s t-test, Figure 1B). Orthologs 


categorized as nonpermissive have levels of infection that are not significantly different from 
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when no DPP4 is present (Figure 1B). Our previous work found that mDPP4 could support 
MERS-CoV infection with just two amino acid substitutions (13) and that knocking out a 
glycosylation site present on blade [TV of mDPP4 was crucial for this result (19). Here, we 
investigate alternate nonpermissive DPP4 orthologs to determine whether they can act as 
backbones to support MERS-CoV infection and if so, what determinants are responsible for 
blocking infection in the wildtype molecules. 

To assess the potential for fDPP4, haDPP4, and gpDPP4 to act as functional receptors for 
MERS-CoV, we compared the overall structures of each molecule. As the crystal structures for 
these orthologs have not yet been solved, we generated homology models using I-TASSER (30). 
Each structure is predicted to have a highly similar backbone topology to hDPP4 (Figure 2A). 
The root mean square deviation (RMSD) scores obtained for mDPP4, fDPP4, haDPP4 and 
gpDPP4 aligned to hDPP4 (PDB code 4L72) are 0.644, 0.616, 0.378, and 0.604, respectively, 
with a smaller number indicating greater structural homology. This can be compared to amino 
acid sequence identity values of 85%, 88%, 85%, and 87%, respectively. The high sequence and 
predicted structural similarities between hDPP4 and these orthologs suggest that they can likely 
act as backbones to support MERS-CoV infection, consistent with previous DPP4 ortholog work 
(15, 17). 

To gain an intuition on specific differences present between permissive and 
nonpermissive DPP4 orthologs, we can assess their relatedness through phylogenetics. 
Unfortunately, analysis of the DPP4 gene tree does not reveal a pattern of permissivity on a 
global sequence scale (i.e., permissive and nonpermissive species do not form distinct 
monophyletic groups) (31-32). However, it is possible that a signature may exist at a local scale. 


Our previous work in mDPP4 revealed that glycosylation on blade IV of the DPP4 molecule can 
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act as a barrier to MERS-CoV infection (19). Thus, we investigated whether glycosylation 
(designated by the motifs NXS or NXT) acts as a broader determinant of DPP4 ortholog 
permissivity. A total of eight N-linked glycans are known to be present in the extracellular 
domain of hDPP4 (29). Without solved crystal structures for DPP4 orthologs, it is unknown how 
many glycans are truly present on the surface of these molecules. However, predictive software 
(33) estimates 7-11 putative glycosylation sites for each of the DPP4 orthologs analyzed in this 
study (Table 1). Aligning the sequences of these orthologs species revealed the presence of 
putative glycosylation sites on blade IV of fDPP4, haDPP4, and gpDPP4 in a region that is 
crucial for interactions with the MERS-CoV RBD (Figure 2B, Table 1). Interestingly, haDPP4 
encodes a glycosylation site at the residues aligning to the site identified in our previous mDPP4 
studies. In contrast, fDPP4 encodes a glycosylation site slightly upstream whereas gpDPP4 
encodes a glycosylation site slightly downstream from the mDPP4 site (Figure 2B). All three 
glycosylation sites lie on blade IV, within a region of DPP4 that interacts with the MERS-CoV 
RBD (Figure 2C). Of note, the downstream gpDPP4 glycosylation site is also present within the 
permissive bDPP4 receptor used in this study (from Pipistrellus pipistrellus), causing us to 
hypothesize that host restriction might be mediated by differences in the sequence or structure of 
gpDPP4 that are independent of this glycosylation site. 

To investigate whether the putative glycosylation sites identified in fDPP4, haDPP4, and 
gpDPP4 play a role in hindering MERS-CoV infection, we assessed receptor permissivity for 
glycosylation knockout mutants. Each DPP4 ortholog knockout includes a mutation that changes 
the N of the glycosylation NXT (or NXS) motif to an alanine, designated by “-gly” (specific 
mutations are f[DPP4 N331A, haDPP4 N332A, and gpDPP4 N548A with residue numbering 


relative to its own sequence). Results show that removing glycosylation from these three DPP4 
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orthologs, confirmed by a ~2.5 kDa downward shift in protein mobility via Western blot analysis 
(Figure 3B), did not result in an increase in infection (Figure 3A). In fact, levels of infection 
supported by all three glycosylation knockout molecules were not significantly different from 
their respective wildtype DPP4 molecules (p < 0.05, Student’s t-test, Figure 3C). 
Immunofluorescence assay (IFA) and flow cytometry confirmed that the DPP4 variants were 
expressed at the surface of the cell (Figure 4), ruling out the possibility of nonpermissivity 
resulting from a lack of proper expression of the DPP4 ortholog glycosylation knockouts. 
Although removing the blade IV glycosylation site on its own is not enough to confer 
permissivity to fDPP4, haDPP4, and gpDPP4 (Figure 3A), we have previously shown in mDPP4 
that a second determinant acts in conjunction with the blade IV glycosylation site to impact the 
interactions between DPP4 and the MERS-CoV RBD (13, 19). Specifically, introducing select 
human amino acids into the mDPP4 sequence on blade V (A288L) acts to support MERS-CoV 
infection when combined with the blade IV glycosylation knockout (N328A). In the case of 
haDPP4, the L288 amino acid identity is conserved between the hamster and human sequences, 
suggesting that different determinants are responsible for blocking MERS-CoV infection 
between haDPP4 and mDPP4. This is consistent with prior studies which identified several 
amino acid differences between hDPP4 and haDPP4 in blade V that were required to support 
MERS-CoV infection (17). Therefore, we asked whether specific changes in blade V of haDPP4 
could allow it to act as a functional MERS-CoV receptor when combined with the blade IV 
glycosylation knockout mutation. We found that three mutations on blade V (E289A, V291M, 
and T293], residue numbering relative to haDPP4) conferred MERS-CoV receptor function to 
haDPP4 when introduced in conjunction with the blade IV glycosylation knockout mutation 


(Figure 5A). Two of the identified mutations overlap with previous work while a third site 


FILSVOMGN JO AINN AQ ZL02 ‘82 Ainr uo /BiowseiAl//:dyjy wo. pepeojumoq 


(0) 
= 
6 
oO 
(0) 
— 
(72) 
O 
jae 
— 
(OL 
= 
O 
72) 
=) 
(= 
= 
TO 
ve 
[ok 
0) 
O 
O 
~s 


Journal of Virology 


Journal of Virology 


191 


192 


193 


194 


195 


196 


197 


198 


199 


200 


201 


202 


203 


204 


205 


206 


207 


208 


209 


210 


211 


212 


213 


knocks out the same blade IV glycosylation site but by mutating a different residue of the 
glycosylation site motif (Table 2). Testing combinations of these mutations in sets that included 
less than the four identified here, including testing haDPP4 E289A, V291M, and T293] (denoted 
haDPP4 (289)) in a glycosylation positive background, did not result in an increase in infection 
(data not shown). These results demonstrate that the presence of the blade IV glycosylation site 
plays an important role in regulating the MERS-CoV receptor function of haDPP4. However, 
species-specific differences in blade V of haDPP4, which differ from those found in mDPP4, are 
also responsible for haDPP4’s inability to act as a functional receptor for MERS-CoV infection. 
Given our findings with haDPP4, we then went on to assess whether the same principles 
applied to f(DPP4. Specifically, we tested whether additional changes in blade V would act in 
concert with the blade IV glycosylation knockout mutation to enhance the MERS-CoV receptor 
function. Previous studies found that swapping fDPP4 amino acids 246 to 503 with hDPP4 
amino acids 247 to 504 allowed fDPP4 to support MERS-CoV infection, but no single mutations 
were identified that could recapitulate the susceptible phenotype (15). Because the mDPP4 and 
haDPP4 data suggest that changes on both blades IV and V of DPP4 are likely required to confer 
infection, we mutated residues on blades V and tried them in combination with the blade IV 
glycosylation knockout mutation (N331A, residue numbering relative to f(DPP4). We found that 
no set of mutations on blade V conferred permissiveness to MERS-CoV when made only with 
the blade IV glycosylation knockout mutant (data not shown). Given this results, we expanded 
our search and made additional mutations on blade IV of fDPP4. By identifying residues not 
present in permissive DPP4 orthologs, we generated a chimeric fDPP4 that changed 10 residues 
on blade IV, including the glycosylation knockout mutation (N331A), and 9 residues on blade V 


to the equivalent human amino acid identities, indicated by their starting residues of 330 and 
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278, respectively, numbering relative to the fDPP4 sequence (Figure 5B). Results show that the 
combination of these 19 amino acid changes on blades IV and V showed a significant increase in 
infection (p < 0.05, Student’s t-test, Figures 5B and 5D). Western blot analysis confirmed high 
expression of each DPP4 variant. The larger downward shift seen for the fDPP4 (278) receptor is 
likely due to the large number of amino acid changes made on blade V (Figure 5C). 
Combinations of individual mutations from these blocks did not result in an increase in infection 
(data not shown). These data reveal that unlike mDPP4 and haDPP4, additional determinants in 
blade IV contribute to the ability of fDPP4 to act as a functional MERS-CoV receptor. 

In addition to haDPP4 and fDPP4, we also investigated the determinants of gpDPP4. 
Whereas haDPP4 and fDPP4 do not share their glycosylation sites with permissive molecules, 
gpDPP4 shares its downstream glycosylation site with bDPP4 (Figure 2B). Because of this, we 
knocked out the blade IV glycosylation site in bDPP4 to determine 1) whether it is truly 
glycosylated and 2) its impact on the permissivity of bDPP4. Results show that removing 
glycosylation from bDPP4 caused no statistically significant change in its ability to support 
MERS-CoV infection (Figures 6A and 6D). Western blot analysis confirmed that the bDPP4 
glycosylation motif is a true glycosylation site when expressed in human cells, as evidenced by a 
downward shift in the glycosylation knockout protein band (Figure 6B). Additionally, flow 
cytometry and IFA confirmed surface expression of each variant tested in this assay (Figures 4B- 
C, 6C). Low apparent expression of bDPP4 is likely due to inefficient binding of the a-hDPP4 
antibody and bDPP4, an observation that is supported by previous work (35). Unfortunately, the 
gpDPP4 plasmid was not stable, and the gene was repeatedly lost during site-directed 
mutagenesis. This instability prevented further investigation of the determinants responsible for 


gpDPP4’s inability to support MERS-CoV infection. Therefore, future studies are needed to 
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determine whether changes on both blades IV and V or just blade IV of gpDPP4 can confer 
permissivity when combined with the blade IV glycosylation knockout mutation, as was shown 


with fDPP4 and haDPP4, respectively. 


DISCUSSION 

The interactions between the virus spike and host cell receptor are crucial for mediating 
infection and for acting as a barrier to nonpermissive species. However, viruses frequently 
overcome species barriers to expand their host range. Glycosylation plays an important role in 
mediating permissivity for many viruses, both in the context of the virus spike protein and of the 
host cell receptor. For example, glycosylation of the HIN1 influenza virus hemagglutinin protein 
is highly conserved and has been shown to be crucial for host cell receptor binding specificity 
(36-37). Conversely, glycosylation of the parvovirus host cell receptor (transferrin receptor type- 
1) in canines was found to confer resistance to infection by previously circulating feline 
parvovirus strains (38). In fact, dogs did not become susceptible to parvovirus until a lineage 
(now canine parvovirus) was able to overcome this glycosylation through specific mutations in 
the capsid protein (38-40). 

Glycosylation has previously been identified as a key determinant of coronavirus host 
range. For example, some lineages of group | coronaviruses use aminopeptidase N (APN) as 
their host cell receptor and APN orthologs have different glycosylation profiles that influence the 
species-specificity of coronavirus infection. Porcine and feline APN orthologs are glycosylated 
and can support infection by porcine and feline coronaviruses, respectively (41). However, 
adding a glycosylation site into human APN near residue 290 can abrogate its ability to support 


human coronavirus 229E (HCoV-229E) infection (41). Similarly, removing the aligning 
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glycosylation site in mouse APN can allow it to act as a functional receptor for HCoV-229E 
(42). As another example, SARS-CoV utilizes angiotensin-converting enzyme 2 (ACE2) as its 
functional receptor. Whereas rat ACE2 is not permissive to SARS-CoV, introduction of a 
glycosylation site combined with a point mutation allows rat ACE2 to support SARS-CoV 
infection (43). Taken together, these results highlight glycosylation as a key determinant of 
receptor species-specificity for numerous coronaviruses. 

Here we investigate the extent to which glycosylation of DPP4 can act to mediate the 
host range of MERS-CoV. Our previous work revealed the importance of glycosylation in 
blocking MERS-CoV infection in the context of mDPP4 (19). Here, we find that glycosylation is 
not the only determinant that mediates MERS-CoV infection in the nonpermissive orthologs 
fDPP4, haDPP4, and gpDPP4. This is particularly surprising in the case of haDPP4; this protein 
not only has a glycosylation site in the same location as mDPP4 (Figure 2B), but has the same 
amino acid identity as hDPP4 at the secondary blade V residue (288) that was identified to be 
important for mDPP4. Instead, we found that three amino acid substitutions in blade V of 
haDPP4 affect the molecule’s ability to mediate MERS-CoV infection when coupled with the 
blade IV glycosylation knockout mutation (Figure 5A). For fDPP4, we found that extensive 
changes on both blades IV and V were required to support MERS-CoV infection (Figure 5B), 
indicating that there are likely species-specific determinants that are not shared with mDPP4 and 
haDPP4. Finally, the blade IV glycosylation site in gpDPP4 is conserved with the permissive 
bDPP4 ortholog (Figure 2B). We found that removing this glycosylation site in bDPP4 did not 
affect its ability to support MERS-CoV infection (Figure 6A). This suggests that the impact of 
the bDPP4 blade IV glycosylation site on receptor permissivity may be species-specific and 


further research should work to elucidate additional changes that allow gpDPP4 to act as a 
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functional receptor for MERS-CoV. These changes can then be made independently and in 
combination with the blade IV glycosylation site to determine its impact within the gpDPP4 
backbone. Overall, we find that glycosylation in DPP4 orthologs is a substantial barrier to 
MERS-CoV infection, particularly when combined with species-specific changes in blades V 
and/or IV. The biochemical mechanism of these additional determinants can be investigated in 
future studies. 


Comparing the data we have so far on the residues that are important for mediating 


permissivity in hDPP4, mDPP4, haDPP4, and fDPP4, it is difficult to discern an obvious pattern. 


The most evident trend is that at least one change is required on both blades IV and V of DPP4 
(Table 2). This indicates that there are two key points of interaction between DPP4 and the 
MERS-CoV RBD which is consistent with previous structural work examining the binding 
interface of MERS-CoV RBD and hDPP4 (29, 34, 44). We show that these two points of 
interaction are also important for allowing the virus to utilize a new species receptor. Further 
understanding of the host range expansion of MERS-CoV will come with identifying which 
mutations can confer permissivity to DPP4 receptors of currently nonpermissive species. 
Additionally, solving the crystal structures of nonpermissive DPP4 orthologs will reveal the 
specific interactions between these molecules and MERS-CoV and help to elucidate what 
prevents MERS-CoV from successfully utilizing these molecules as functional receptors. The 
inability to identify a small number of changes that confer permissivity to DPP4 orthologs 
suggests that extensive remodeling of the MERS-CoV RBD might be required for the virus to 
infect these nonpermissive species. Access to the crystal structures of these DPP4 orthologs 
would enable us to model changes in the MERS-CoV RBD that would allow the virus to 


successfully infect nonpermissive species, thereby enhancing our ability to study MERS-CoV 
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pathogenesis and provide robust small animal models for evaluating vaccines and therapeutics. 
As of now, our results indicate that the generation of a transgenic small animal model in these 
alternate species would require extensive genomic editing due to the number of mutations 
needed for each DPP4 to support MERS-CoV infection. 

The importance of glycosylation in blocking MERS-CoV infection may vary between 
species. To gain a better intuition about the extent of glycosylation among DPP4 orthologs, we 
constructed a phylogenetic tree of a subset of full-length DPP4 protein sequences (Figure 7). In 
the tree, shaded colors indicate the general organism group that each species belongs to - blue: 
reptiles and amphibians; green: avian species; orange: other mammals; red: Chiroptera (bats); 
purple: ungulates; gray: rodents; pink: primates (Figure 7). The DPP4 gene tree is slightly 
discordant with the species tree, notably with the horse and African savanna elephant DPP4 
sequences not clustering with other ungulate (purple) orthologs. Plotted adjacent to the 
phylogenetic tree are glycosylation sites that are either upstream (column 1), at the same site 
(column 2), or downstream (column 3) of the glycosylation site that is present in mDPP4 
(residues 328-330). Permissivity data is indicated in the far right column, with green squares 
indicating permissive species and red squares indicating nonpermissive species, based on either 
in vitro or in vivo data (10-17). 

The small number of known permissive and nonpermissive species makes it difficult to 
map potential shifts in permissivity on the phylogenetic tree. Because permissive and 
nonpermissive species do not cluster together on the DPP4 tree (Figure 7), it is more likely that 
local differences between each protein is more important than global similarity. Despite our 
inability to identify a distinct pattern, a few interesting observations do emerge. First, all non- 


human primates lack glycosylation sites near residue 330 (equivalent to residue 336 in humans) 
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(Figure 7). Current data would postulate that all non-human primates are susceptible to MERS- 
CoV infection, but testing this hypothesis in vitro would help reveal whether any of these 
orthologs are nonpermissive. If so, a mechanism other than glycosylation at our location of 
interest would be responsible for blocking MERS-CoV infection. Second, whereas other 
permissive DPP4 orthologs have glycosylation sites in this region (horse and common pipistrelle 
DPP4), of note is the lack of a glycosylation site in the nonpermissive domestic pig, sheep, and 
cattle DPP4 molecules (Figure 7). This suggests that glycosylation in this region is not the 
primary explanation for why these species do not support MERS-CoV infection. One example 
could be the incompatibility of host cell proteases which are required for proper cleavage of the 
MERS-CoV spike prior to cell entry. Previous works has shown that spike cleavage status 
determines cell tropism (45) and that adaptation to cellular proteases may have played an 
important role in the cross-species transmission of MERS-CoV from bats to humans (35). 
Follow-up studies could investigate the mechanism of nonpermissivity in these species 
specifically and whether or not it lies at the level of the receptor. Third, glycosylation sites in the 
designated region are prevalent in rodents, other mammals, and avian species (Figure 7). For 
example, eight other species have putative glycosylation sites aligning to the glycosylation site 
identified in mDPP4. The majority of these are present in the rodent (gray) group, however two 
are present within the Chiroptera group (black flying fox and large flying fox) (Figure 7). 
However, the diversity of glycosylation profiles suggests a lack of strong conservation of the site 
seen in mDPP4. Further research to determine whether these other orthologs are permissive can 
help elucidate whether this region plays a broader role in the MERS-CoV infection phenotype. In 
general, while some trends are seen when looking at the glycosylation profile across many 


species (e.g. upstream glycosylation sites in the other mammal (orange) group), more data on 
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species permissivity will help determine whether DPP4 phylogenetic relationships can help 
inform receptor binding dynamics. 

Here, we demonstrate that glycosylation is an important barrier to MERS-CoV infection, 
yet other species-specific determinants are also responsible. Removing the glycosylation sites of 
fDPP4, haDPP4, and gpDPP4 orthologs did not result in an increase in infection. Rather, a set of 
19 and 4 amino acid substitutions in fDPP4 and haDPP4, respectively, were required before a 
significant increase in infection could be detected. Future work can focus on revealing the 
importance of these additional substitutions in mediating MERS-CoV infection. By further 
elucidating the mechanisms by which DPP4 orthologs block MERS-CoV infection, we can better 
understand the constraints on coronavirus host range expansion and predict coronavirus 
emergence into new species in the future. Furthermore, this work can inform the development of 
better animal models for studying MERS-CoV pathogenesis. Our previous work on mDPP4 
helped generate a successful transgenic mouse using the CRISPR-Cas9 gene editing technique 
(24). The success of this model emphasizes the importance of understanding receptor-virus 


interactions and is a potent example of the direct application of these studies. 


MATERIALS AND METHODS 
Viruses and cells 

MERS-CoV was isolated from a molecular clone as described in Scobey et al. (46). 
Recombinant MERS-CoV tagged with tomato red fluorescent protein ((MERS-CoV-RFP) was 
engineered based on the EMC2012 substrain and was shown to infect and replicate similarly 
compared to wildtype MERS-CoV (46). Recombinant viruses were passaged once (P1) on Vero 


cells to generate a working stock of (MERS-CoV-RFP. All recombinant viruses were produced 
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and studied under BSL3 conditions, wearing Tyvek suits, apron and booties, PAPR, and double 
gloved as described in the Laboratory Safety Plan. 

Measures of rMERS-CoV-RFP infectivity, receptor protein expression levels, and 
immunofluorescence assays (IFA) were performed in human embryonic kidney (HEK) 293T 
cells transfected with cDNAs encoding wildtype or mutant DPP4 orthologs. Cells were seeded at 
1 x 10° cells per well in 6-well plates pretreated with 1 mL of 10 pg/mL poly-L-lysine (Sigma) 


in 1X DPBS (Gibco). All cells were grown at 37 °C with 5% COd. 


DPP4 ortholog constructs 

DPP4 orthologs were ordered on plasmids (cDPP4, hDPP4, mDPP4) or constructed using 
Gibson assembly (bDPP4, fDPP4, haDPP4, gpDPP4) (47). DPP4 orthologs were transferred into 
the 945ARRE expression vector, a lentiviral vector derived from pTK945, using the NotI and 
Spel restriction sites. DPP4 mutants were generated by overlap PCR mutagenesis and verified by 
Sanger sequencing. All plasmids were stored at -20 °C. Accession numbers for DPP4 sequences 
used in this study are as follows: bDPP4, KC249974.1; cDPP4, XM_006176808.2; fDPP4, 
DQ266376.1; gpDPP4, XM_013142395.1; haDPP4, NM_001310571.1; hDPP4, NM_001935.3; 


mDPP4, NM_010074.3. 


DPP4 ortholog transfections 

DPP4 orthologs were transfected into HEK 293T cells, as these cells express little if any 
endogenous hDPP4 (13). Transfections followed a 2X BES protocol. For each DPP4 ortholog, 
15 pg of the construct was mixed with 15 ug of pecDNA3.1, 125 uL of 1 M CaClo, and dH20 up 


to final volume of 500 uL. Then, 500 uL of 2X BES was added to each sample in a dropwise 
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manner, vortexing between drops. Samples were incubated at room temperature for 45 min. 
During this incubation, media was removed from 6-well plates seeded with 1 x 10° HEK 293T 
cells/well and 1 mL of fetal bovine serum (FBS)-free media was added to each well. After the 45 
min incubation, 200 uL of each transfection reaction was added to individual wells, dropwise 
and in concentric circles. Plates were incubated at 37 °C for 2 hours before 150 uL of FBS was 
added to each well and the cultures incubated overnight at 37 °C. To measure transfection 
efficiencies, venus protein-tagged 945 plasmid was transfected in parallel with each experiment 


(with efficiencies typically between ~80-100%). 


Infection of DPP4-transfected cells by MERS-CoV 

At ~20 hours post-transfection, media was removed and each well infected with 200 uL 
of rMERS-CoV-RFP at a multiplicity of infection (MOD of 5 for 45 min at 37 °C. Next, 2 mL of 
3% FBS media was added to each well and the plates incubated at 37 °C for 24 hours. Cells were 
visualized 24 hours post-infection in the BSL3 using an Olympus IX51 microscope at a 
magnification of 10X and wavelength of 541 nm. Each well was imaged using a Hamamatsu 
ORCA-Flash4.0 LT camera. Successful virus infection is indicated by the presence of red 


fluorescence. 


Fluorescent cell count assay 

For quantification of levels of MERS-CoV infection, transfected cells were infected as 
above, except at an MOI of 0.1. After the 45 min incubation period, virus samples were removed 
and the cells washed twice with | mL 1X DPBS. Next, 2 mL of 3% FBS media was added to 


each well and infections incubated at 37 °C for 72 hours. At 72 hours post-infection, the number 
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of red fluorescent cells was determined. Briefly, three randomly selected microscope fields were 
counted for each well. The average number of red cells was taken from these three fields and 
multiplied by a factor of 537 to approximate the number of red cells in the entire well. Infection 
in the presence of each DPP4 ortholog was measured as described in triplicate, thus 9 fields were 


measured to determine significance between two DPP4 treatments using a Student’s t-test. 


Protein analysis and Western blotting 

Cells were transfected with DPP4 constructs as described above. At 48 hours post- 
transfection, media was removed from each well and cells were washed once with 2 mL 1X 
DPBS. Next, 200 uL of ice cold AV lysis buffer (20 mM Tris-HCl, 150 mM NaCl, 0.5% 
deoxycholine, 1% Igepal, 0.1% SDS) was added to each well. Cells were incubated at room 
temperature for 10 min with the lysis buffer. Lysates were then transferred to 1.5 mL 
microcentrifuge tubes and centrifuged at maximum speed for 2 min. Supernatants were 
transferred to new tubes and 200 pL of a 1OmM EDTA/0.9% SDS stock added to each tube. 
Lysates were stored at -80 °C. 

Protein concentration for each sample was measured using a Bradford assay. Briefly, an 
aliquot of each sample was diluted 1:5 and 10 uL added to a 96-well plate in triplicate. Bovine 
serum albumin (BSA) was used as a standard with 10 uL of a dilution series added to the 96-well 
plate in triplicate. To each well, 200 uL of Coomassie Plus protein assay reagent 
(ThermoScientific) was added and the plate incubated at room temperature in the dark for 10 
min. Absorbance readings were measured at a wavelength of 595 nm and the protein 


concentration of each sample determined by comparison to the BSA dilution standard curve. 
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Proteins were separated using either a 6% (for DPP4) or an 8% (for actin) polyacrylamide 
gel. For each sample, 15 ug of the protein lysate was mixed with 6X loading dye and heated to 
95 °C for 10 min. Samples were loaded into the gel and separated by electrophoresis in 1X 
running buffer (3.03g Tris HCI, 14.4g glycine, 1g SDS in 1L dH2O) at 55 mA for ~2 hours. 
Proteins were transferred onto an Immun-Blot PVDF membrane (BioRad) at 15V for 45 min 
using 1X dry transfer buffer (5.82g Tris base, 2.92g glycine, 100 mL methanol in 1L dH20). The 
membranes were blocked in 1X phosphate-buffered saline with Tween 20 (PBST) with 5% milk 
at room temperature for 1 hour with shaking. Membranes were then treated with 1° antibody: 
goat a-DPP4 (R&D systems) at 1:1000 or goat a-actin (Santa Cruz Biotechnology) at 1:1000 in 
1X PBST and incubated at 4 °C overnight with shaking. Membranes were washed three times 
with 1X PBST (10 min per wash) before treatment with 2° antibody. Both DPP4 and actin blots 
used rabbit a-goat HRP at 1:10,000 in 1X PBST with 5% milk and were incubated at room 
temperature for | hour with shaking. Membranes were washed three times with 1X PBST (10 
min per wash) and developed using Amersham ECL Western Blotting detection reagents (GE 
Healthcare). Membranes were incubated at room temperature for 10 min before imaging using 


Amersham Hyperfilm ECL (GE Healthcare). 


Immunofluorescence Assay (IFA) 

For IFA analyses, transfections were performed in 50 mm glass-bottom dishes seeded 
with 2 x 10° HEK 293T cells. At 24 hours post-transfection, cells were washed once with 1X 
DPBS and then fixed with 2 mL 2% paraformaldehyde (PFA) in 1X DPBS for 30 min at room 
temperature. Cells were washed once with 2 mL of 100 mM glycine in 1X DPBS and incubated 


in 2 mL of 100 mM glycine in 1X DPBS for 15 min at room temperature. Cells were washed 
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twice with 1X DPBS before blocking. To each plate, 200 uL of 3% BSA and 10% normal 
donkey sera in IFA wash (0.05% Tween-20 in 1X DPBS) was added and incubated at room 
temperature for | hour with shaking. Next, cells were incubated with goat a-DPP4 polyclonal 
antibody (R&D systems) at 1:50 in 3% BSA in IFA wash at room temperature for | hour with 
shaking. Cells were washed three times with IFA wash and then incubated with 2° donkey-anti- 
goat Alexa Fluor 488 (Life Technologies) at 1:500 for 1 hour at room temperature in the dark 
with shaking. Cells were washed again three times with IFA wash. To each plate, three drops of 
ProLong Gold antifade reagent (Invitrogen) were added to the cells before imaging at 20X 
magnification. Green fluorescence indicates DPP4 staining while blue is DAPI staining for 


nuclei. 


Flow cytometry analysis 

Transfections were performed as detailed above. At ~20 hours post-transfection, media 
was removed and cells were washed with 1X DPBS. Cells were removed from plates using non- 
enzymatic cell dissociation buffer (Gibco). Cells were then resuspended in complete media and 
centrifuged (all spins at 500 x g for 5 minutes). Supernatants were removed and cells were 
washed three times with 1X DPBS. Cells were counted and resuspended at a concentration of 5 
x 10° cells/mL. To each well of a 96-well plate (Corning), 200 uL of cell suspension was added 
in duplicate, plates were spun to pellet cells, and the supernantant was discarded. Cells were 
suspended in 100 uL Flow Cytometry Staining Buffer (eBiosciences) with 8 ug/mL goat a-DPP4 
polyclonal antibody (R&D Systems) and incubated at 4 °C for 1 hour. Cells were washed twice 
with Staining Buffer and then incubated with | uL/well of FITC-conjugated donkey a-goat 


secondary antibody (Life Technologies) in 100 pl Staining Buffer for 1 hour at 4 °C. Cells were 
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washed twice with Staining Buffer, and resuspended in 150 uL Staining Buffer. At least 10,000 
events per well were acquired on the Guava easyCyte HT flow cytometer (Millipore) for each 


read. The background signal of secondary antibody only (no primary antibody) was used to 


subtract out the background. Data were analyzed using FlowJo software version 10.3 (Tree Star). 


Phylogenetic analysis 

Amino acid sequences were retrieved from GenBank and aligned using MAFFT (48). 
The phylogenetic tree was generated using maximum likelihood and 100 bootstrap replicates 
with the PhyML package. Smart model selection (SMS) was implemented to determine the best 
model parameters (49). The tree was visualized using Evol View (50). Only bootstrap values >50 


are displayed. 


Statistics 
All quantitative data are presented as means + one standard deviation. All statistics were 


performed using the R programming language version 3.1.2. 
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FIGURE LEGENDS 


FIG 1. Permissivity of DPP4 orthologs to MERS-CoV. (A) Seven DPP4 orthologs were tested 
for their ability to support infection by r'MERS-CoV-RFP. DPP4 constructs were transfected into 
HEK 293T cells and infected at an MOI of 5 ~20 hours post-transfection. Cells were imaged for 
fluorescence ~24 hpi. hDPP4, human DPP4; cDPP4, camel DPP4; bDPP4, bat DPP4; mDPP4, 
mouse DPP4; fDPP4, ferret DPP4; haDPP4, hamster DPP4; gpDPP4, guinea pig DPP4. (B) 
Mean fluorescent cell counts of MERS-CoV infection utilizing various DPP4 orthologs. Cells 
were infected at an MOI of 0.1 and the number of infected cells counted at ~72 hpi. Each DPP4 
ortholog was measured in triplicate. Only hDPP4, bDPP4, and cDPP4 have significantly higher 


levels of infection compared to no DPP4 (p < 0.05, Student’s t-test, indicated by *). All DPP4 
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orthologs have significantly lower levels of infection compared to hDPP4 (p < 0.05, Student’s t- 
test). Levels of infection are not significantly different between bDPP4 and cDPP4. Error bars 


indicate mean values + one standard deviation. 


FIG 2. Sequence and structural comparison of nonpermissive DPP4 orthologs. (A) Structural 
comparison of threaded molecules (30) mDPP4 (orange), fDPP4 (green), haDPP4 (blue), and 
gpDPP4 (purple) overlaid on hDPP4 (yellow) complexed with the MERS-CoV RBD (red) (PDB 
code 4L72). (B) Sequence alignment of permissive (human, camel, bat; blue) and nonpermissive 
(mouse, ferret, hamster, guinea pig; red) DPP4 amino acid sequences. Residue 330 is numbered 
relative to mDPP4. Boxes represent glycosylation sites that are either unique to nonpermissive 
species (black) or shared with a permissive species (gray). (C) hDPP4 (yellow) complexed with 
the MERS-CoV RBD (red) (PDB code 4L72). Residues aligning to the ferret (green), hamster 
and mouse (blue), and guinea pig (purple) glycosylation sites are highlighted. Dashed line circles 


indicate the regions of the DPP4 molecule that corresponds to blades IV and V. 


FIG 3. DPP4 ortholog glycosylation knockout mutants. (A) Neither wildtype nor glycosylation 
knockout DPP4 molecules for ferret (f(DPP4), hamster (haDPP4), or guinea pig (gpDPP4) 
support infection by MERS-CoV. (B) Successful removal of glycosylation is supported by a ~2.5 
kDa downward shift seen via Western blot. Top blot represents DPP4 and the bottom blot 
represents B-actin as a control. (C) Fluorescent cell counts of MERS-CoV infection utilizing 
various DPP4 orthologs and their respective glycosylation knockout mutants. Cells were infected 
at an MOI of 0.1 and number of infected cells counted at 72 hpi. Each DPP4 ortholog was 


measured in triplicate. Only hDPP4 has significantly higher levels of infection compared to no 
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DPP4 (p < 0.05, Student’s t-test, indicated by *). The remaining DPP4 orthologs have infection 
levels that are not significantly different from no DPP4 or each other. Error bars indicate mean 


values + one standard deviation. 


FIG 4. DPP4 and mutant variants are expressed on the surface of cells as evidence by 
immunofluorescence assay (A) and flow cytometry (B, C). (A) Cells were transfected with each 
DPP4 ortholog, fixed, and probed with primary goat-anti-DPP4 polyclonal antibody (R&D 
Systems) at 1:50 and secondary donkey-a-goat Alexa Fluor 488 (Life Technologies) at 1:500. 
Cells were imaged at 20X for DAPI (300 ms exposure) and DPP4 (1.5 s exposure). (B) DPP4 
expression frequencies (blue outlined histogram) by DPP4 construct after subtraction of 
background from replicate wells stained with secondary donkey-a-goat IgG (H+L) Alexa Fluor 
488 antibody only (grey shaded histogram). Percentage values represent the average across two 
duplicate wells. (C) Geometric mean fluorescence intensity of the DPP4 positive populations for 


each DPP4 construct. 


FIG 5. Many amino acid changes are required to make fDPP4 and haDPP4 permissive to 
MERS-CoV infection. (A) Removing glycosylation on its own does not confer permissivity to 
haDPP4. However, combining three amino acid changes on blade V (starting at residue 289) 
with the glycosylation knockout mutant on blade [TV (N332A) results in high levels of MERS- 
CoV infection. Sequences show the alignment between hDPP4 and haDPP4 with the blue boxes 
indicating the amino acids that were swapped from hDPP4 into haDPP4. (B) Removing 
glycosylation on its own does not confer permissivity to fDPP4. However, introducing a set of 


amino acid changes on blade V (starting at residue 278) and blade IV (starting at residue 330) 
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allows fDPP4 to support MERS-CoV infection (f(DPP4 (278)(330)). Sequences show the 
alignment between hDPP4 and fDPP4 with the blue boxes indicating the amino acids that were 
swapped from hDPP4 into fDPP4. Note that fDPP4 —gly is a negative control and only includes 
the single point mutation N331A. (C) Western blot analysis of fDPP4 and haDPP4 and 
designated variants for DPP4 and B-actin expression. Successful glycosylation knockout is 
indicated by a downward shift of ~2.5 kDa. (D) Fluorescent cell counts of MERS-CoV infection 
utilizing DPP4 orthologs. Cells were infected at an MOI of 0.1 and number of red cells counted 
at 72 hpi. Each DPP4 ortholog was measured in triplicate. hDPP4, fDPP4 (278)(330), and 
haDPP4 (289),-gly have levels of infection that are significantly greater than no DPP4 (p < 0.05, 
Student’s t-test, indicated by *). fDPP4 —gly and haDPP4 —gly infection levels are not 


significantly different from no DPP4. Error bars indicate mean values + one standard deviation. 


FIG 6. Bat and guinea pig DPP4 share the same glycosylation site downstream of the site 
identified to be important in mDPP4 (Figure 2B). (A) Removing the glycosylation site from 
bDPP4 shows no decrease in infection while removing glycosylation from gpDPP4 resulted in 
no increase in infection. (B) Western blot analysis of bDPP4 and gpDPP4 and their respective 
glycosylation knockout mutants for DPP4 and B-actin expression. Successful glycosylation 
knockout is indicated by a downward shift of ~2.5 kDa. (C) Fluorescent cell counts of MERS- 
CoV infection utilizing various DPP4 orthologs. Cells were infected at an MOI of 0.1 and 
number of infected cells counted at 72 hpi. Each DPP4 ortholog was measured in triplicate. 
hDPP4, bDPP4, and bDPP4 — gly have levels of infection that are significantly higher than no 
DPP4 (p < 0.05, Student’s t-test, indicated by *). goDPP4 and gpDPP4 —gly infection levels are 


not significantly different from no DPP4. (D) DPP4 and mutant variants are expressed on the 
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759 surface of cells, visible by immunofluorescence. Cells were transfected with each DPP4 

760 ortholog, fixed, and probed with primary goat-anti-DPP4 polyclonal antibody (R&D Systems) at 
761 ~—‘1:50 and secondary donkey-anti-goat Alexa Fluor 488 (Life Technologies) at 1:500. Cells were 
762 imaged at 20X for DAPI (300 ms exposure) and DPP4 (1.5 s exposure). Error bars indicate mean 


763 values + one standard deviation. 
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765  FIG7. DPP4 protein phylogenetic tree based on amino acid sequences. Shaded colors indicate 
766 the group each species falls in. Blue: reptiles and amphibians; Green: avian species; Orange: 
767 other mammals; Red: Chiroptera (bats); Purple: ungulates; Gray: rodents; Pink: primates. 

768 Colored circles to the right of the species names indicate whether the sequence has a 

769 glycosylation site upstream (first column), at the same site (second column), or downstream 
770 (third column) of the NXT glycosylation site in mDPP4 (residues 332-334). Numbers inside the 


771 circle designate how many amino acids upstream (or downstream) the N of the NXT or NXS 


Journal of Virology 


772 glycosylation site is. For the second column, a circle indicates that there is a glycosylation site 
773 aligning to the site present in mDPP4. Squares in the rightmost column indicate permissive 


774 (green) or non-permissive (red) species, as determined from either in vivo or in vitro studies. 
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hDPP4 85 | 92 | 150 219 | 229 281 | 321 520 
bDPP4 82 | 89 | 147 216 | 226 | 272 315 332 | 393 | 432 | 490 
mDPP4 83 | 90 | 144 223 315 | 328 514 
fDPP4 84 | 91 178 228 280 | 320 437 519 
haDPP4 83 | 90 | 148 227 275 319 | 332 518 
gpDPP4 | 32 | 295 | 302 | 360 439 491 | 531 548 730 

777 


778 Table 1: Known and putative N-linked glycans in DPP4 orthologs 
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Numbers indicate the residues known (hDPP4) or predicted (all other DPP4 orthologs) (33) to be 
N-linked glycosylated for each DPP4 molecule. All residues are numbered relative to the given 
DPP4 molecule, with columns representing residues that align. Out of the eight N-linked glycans 
in hDPP4, only four are conserved across all orthologs. Bold font indicates sites that align to 


within two residues of a contact residue in hDPP4 (see 29). 


Ortholog Important residues Source 
Human 
267 wns | 295) 336 4 
DPP4 a 
Mouse 
294 336* 13, 24 
DPP4 ere) 
Hamster 291 295 336* | 341 | 346 (16, 17) 
DPP4 293 295 || 297 ean This study 
Ferret 279 - 295 331-341 This study 
DPP4 246 - 505 (15) 


Table 2: Residues identified to be important for MERS-CoV permissivity 


Experimental studies and the residues that have been identified as important for mediating 
permissivity to MERS-CoV among various DPP4 orthologs in vitro. All residues are relative to 
the aligning residue in hDPP4. Residues are either on blade IV (blue shading) or on blade V 
(gray shading). * indicates residues that knock out a glycosylation site in the mouse or hamster 


DPP4 molecules. 
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